: arXi V : astro-ph/0 10623 Ov 113 Jun 2001 



Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



The velocity field of collapsing spherical structures 
Limitations of the spherical infall model in mass estimation 

N. Hiotelis 

Lyceum of Naxos, Chora of Naxos, Naxos 84300, Greece 
email: hiotelis@avra.ipta.demokritos.gr 

Received 15 January 2001 / Accepted 28 May 2001 

Abstract. We assume that the amplitude of the caustics in redshift space is a sum of two components: the first 
one can be predicted by the spherical infall model with no random motion, and the second is due to the random 
motion distribution. Smooth model curves are used to estimate the maximum values of the first component for 
the Coma cluster. Then, an approximation of the radial component of the infall velocity -based on the above 
curves- is derived and a mass profile of the cluster is calculated. This mass profile, that is an upper limit for the 
spherical infall model, combined with estimations given by other authors provides an approximation of a lower 
limit for the mass of the system. 
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1. Introduction 

The spherical infall model has been extensively discussed 



in th e literatur e (Gunn fc Gott [1972| ; Silk |1974| ; Gunn 
19771 ; Peebles ( |1976| , |1980D ; Schechter |1980| ; Schectman 



1982; Ostriker et al. 1988). This model assumes that 



galaxy clusters started as small density perturbations in 
the early universe. These perturbations eventually devi- 
ate from the general expansion and after reaching a maxi- 
mum radius they start collapsing. The most probable sce- 
nario is that galaxies formed during the expansion phase 
of their cluster. Thus, the cluster consists of largely indi- 
vidual galaxies during its collapse. 

The observed velocities of these galaxies (along the 
line-of-sight) give important information about the dy- 
namical state of the cluster (Kaiser 1987). Regos & Geller 
( |1989D showed that plotting the observed velocities of 
galaxies as function of their angular distances from the 
centre of the cluster one obtains a velocity distribution 
that is bounded by sharp, characteristic trumpet-shaped 
curves. These curves, called caustics, form an envelope 
containing the galaxies of the cluster. The caustics are 
used to estimate the d ensity parameter of the universe 
(Regos & Geller 1989 ) and t he mass profile for t he ga laxy 
clusters (Geller et al. 199E ; Reisenegger et al. 2000 ). In 
Sect. 2, we show that the form of the caustics fully defines 
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the profile of the radial infall velocity in the case of a pure 
spherical collapse. In Sect. 3, we describe the spherical 
model. This model is applied in Sect. 4 to data from the 
Coma cluster. In Sect. 5, the conclusions are given. 

2. Observed velocities 

For the sake of simplicity, let us consider the xy plane of a 
spherically symmetric cluster, undergoing a radial collapse 
(Fig.|l]), where r is the position vector of a galaxy at K, 
relative to the centre C of the cluster. The velocity v of 
the galaxy at K relative to the observer is given by the 
relation 
dv 

- = ^+-c, (1) 

where vq is the velocity of the centre C relative to O. 
Defining Vr{r) = dr/dt, where r is the magnitude of r 
and multiplying both sides of eqn.(|l|) by the unit vector 
fx along Oft direction, we have 



Vohsir,R) ^ vc ■ rK± \ I 



7 ) Mr), 



(2) 



where Wobs — ■ ^'k is the observed velocity along the 
line-of-sight. From eqn.(^, for r = R and for r = tm, the 
observed velocity takes the same value Uobs('', ^) = ■ 
f"K. Thus, Vohs takes extreme values for r between R and 
tm- We denote by tm the radius of maximum expansion 
(turnaround radius) where Vrirm) = 0. 
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For constant R, the quantity vc ■ does not depend 
on r and thus the condition dvohs/dr = for the extreme 
is written 



1 - 



dvr{r) 
dr 



— «,.(r)^0. 



(3) 



The solution r» for r of eqn.(|^) is a function of R (r* = 
rt{R)). This is the radial -relative to the centre of the 
cluster- distance of the galaxy which shows the observer 
the maximum or minimum velocity at a given R. Thus, 
the extreme values of Vohs at R are 



Vohs,cxtiR) = f c • ■tk ± \ / 1 



R 



r.iR) 



VriMR)). (4) 



The curves described by eqn.(^ are known as caustics 
(Regos & Geller ([T989|)). 



Following Diaferio et al. ( 1997 ), we define the ampli- 
tude A of the caustics by the relation 



MR) = ^ bobs, max (i?) - Wobs.minli?)], 

which is given by 



(5) 



^(i?) = -1 



1 - ( ^ ) Vr{r^). 



(6) 



Note that by definition Vr is negative in the infall region 
of the cluster. The amplitude contains only information 
about the dynamical situation of the cluster, since effects 
due to the motion of the observer relative to the cluster are 
excluded. Such effects are described by the term vc ■ vk 
in eqn. and have been studied by Praton & Schneider 
( |1994D . Differentiating eqn.(||) with respect to R, we have 
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Fig. 1. Geometry of the problem. O: observer, K: galaxy, 
C: cluster's centre 



cL4 



Rvrir^) 



1 - 



R^Vr[r.^ 



(5)- 



R \ dvr (?■* ) 
/ dr* 



dr* 
dR' 



(7) 



Because of eqn. (^, the term inside brackets equals zero 
and this leads to 



vl{r;) _ A(R) dA(R) 
rl ^ R dR ' 

Using eqns. (o) and (ph we have 



ddn_A(i?)\~^ 
d In i? / 



(8) 



(9) 



Solving eqn.(g) for R and substituting in eqn.(g) the radial 
velocity is fully defined. This is an interesting result, since 
it shows that the exact data can give the exact profile 
of the infall velocity, in the case of a completely radial 
collapse. It is interesting to note that the above relations 
lead to eqn.jiol), relating the logarithmic derivatives of the 
amplitude to the infall velocity; 



dln[-Vr{r^)] _ dlnA(i?) 



din r* 



dlni? 



(10) 



The above equations refer to the ideal case of a pure spher- 
ical collapse and completely accurate observations. In real 
systems, the velocity field is a superposition of a radial 
systematic component and a component of a random na- 
ture. The first one can be assumed as spherically sym- 
metric while the second accounts for the effects of small- 
scale substructure and observational errors. The effects of 
small-scale substructure are clearly shown in the results of 
N-body simulations performed by van Haarlem and van de 
Weygaert ( 1993| ). They conclude that the velocity profile, 
as predicted by the spherical infall model, compares badly 
with the actual velocity field resulting from their simula- 
tions. In real systems, such as the Coma cluster, the smear- 
ing of the form of the caustics is also clear (van Haarlem 
et.al 1993). However, if the systematic radial component 
of the velocity field is known, a mass profile of the in- 
fall regions of the clusters can be determined, applying 
the spherical infall model. This profile should give a lower 
limit for the mass of the system since the effects of small- 
scale substructure reported above increase the amplitude 
of the caustics. Thus, we apply the previous eqns. assum- 
ing that Vr describes the systematic radial component of 
the velocity. In Sect. 4 we use model functions for the am- 
plitude A, calculating the profile of by eqns. (^) and (^, 
we apply the spherical infall model and we derive a mass 
profile for the Coma cluster. These model functions for A 
are smooth and decreasing, as it required by eqns. (|[) and 
(^, and their approximation is based on the respective 
curves estimated by Geller et al. (1999). 
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3. The spherical infall model - Mass estimation 



The eqn. of motion of a spherical shell of initial radius a 
is given by Newton's law as 

.. _ d^r _ GM{r) 



(11) 



where M (r) is the mass inside radius r and G is the gravi- 
tational constant. If the mass inside the shell is conserved 
(no shells crossing), then M{r) is constant and equal to 
M{a). We assume that this condition holds in what fol- 
lows. Multiplying eqn.(|ll|) by f and integrating, we have 



v^{t) =2E + 



2GM{a) 



(12) 



where E = 1/2 — {GM{a))/a is a constant with di- 
mensions of specific energy and Vi is the initial velocity. 
The solution of the above eqn., in the case of negative E, 
is given in parametric form by the expressions 



GM{a) 
GM{a) 



(1 — COST]) 



(ry - sinrj) + ta- 



(13) 



(14) 



In order to synchronize all shells, we add the constant ta 



given by ta^U^ 



GM[a 



^ {rji — sin rji), where ti is the age of 



the universe at the initial conditions and T^j is the initial 
phase of the shell, given by eqn.(|l3|) for r = a. A shell 
reaches the radius of maximum expansion for rj — n and 
completes its collapse for i] ^ 2tt. We denote by fi{a) the 
fractional excess of mass of a sphere of radius a, relative to 
the mass Mf,(a) of a sphere of the same radius that has a 
constant density (equal to the mean density of the universe 
at the initial conditions). Using the relation between the 
mean density of the universe pb, the Hubble constant H 
and the density parameter fl, that is 



Pb{t) = 



8^ 



(15) 



(16) 



the mass inside radius a is given by 

M(«) = ^^(l + M(a)), 

where index i stands for the values at the initial condi- 
tions. Assuming that the initial velocity of the shell is 
that of Hubble's flow {vi ~ aHi), the energy E is written 



E = ^H^a^n, [n,-^ - 1 - fi{a) 



(17) 



Substituting eqn.(p7|) in eqns. ([l3|), ([TJ) and calculating 
the time derivative of the radius yields the relations 



afieia)+Qi ^ , , 

o n (1 - COST?) 

2 Me (a) 



(18) 



2Hin^^fiJ{a) 



{t] - sin 77) -I- ta 



Vr = aHiil,'^ fij (a) ^ ^^"^^ (20) 
1 — cos 77 

where p,e{a) = 1 + u{a) — fl^^- 

Combining eqns. (^8|), (19) and (^) results in 

it-taf-^'^p^l^^, (21) 

r (1 — cos?7)^ 

where t is the present age of the shell. The time t is ap- 
proximated hy t = Hq~^ f{Qo, z), where z is the redshift of 
the cluster, Hq is the present value of the Hubble constant 
and f2o is the value of the density parameter. The function 
/ for the Friedmann-Lemaitre cosmological model is given 
in sta ndard books of cosmology (see Zel'dovich & Novikov 
|l983 ) and its form depends on the value of Qq relative to 
the critical value fJo = 1. Reasonable initial conditions 
are those at the epoch of decoupling, when matter starts 
playing a significant role. The value of ta is very small 
compared to t, so it can be omitted without significant 
error in the calculation of the phase r/. However, eqn. (21) 
can also be solved in its full form based on the following 
procedure. Solving eqn. (p^) for /ie gives 



Me (a) 



a(l — cosri)fli 



(22) 



2r — a{l — cos 77) 

The substitution of the above expression in eqn. (H) leads 
to the following cubic for a 

+ {1- cos 7j)^Ca - 2r(l - cos 77)C = (23) 

where C = 7;^/(sin^ vHf)- This has one real and positive 
root, given by the relation 



(24) 



Hf{l + cosri) 



(l + VF)^ + {1-VF) 



where 



F = 1 



1 (1 — COSTff 
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(25) 



Hf{l + con riY 

Thus, Me can be expressed as a function of r,Vr, Hi^fli 
and r/. Then, ta can also be expressed in terms of the 
same variables, since it is given by the relation 



ta — ti 



1 fJ.e{a)+^, 



2HM? 



mI (a) 



-V) 



(26) 



2^1 (a) 



and ti is the age of the universe 



where (f ■ 

(/Je(a)+nr^) 

at the initial conditions (decoupling), given i n stan dard 
books of cosmology (see Zel'dovich & Novikov 1983 ). 

The present method can be summarized as follows: A 
known profile of the amplitude of the caustics, using eqns. 
(||) and (^) leads to the profile of radial velocity. Then, 
the solution of eqn. (pT|) defines the phase 77 of any shell. 
Assuming no shell crossing during the evolution of the 
cluster, the combination of eqns. (|l6|), (18) and (|2^) gives 
the mass inside a spherical region with current radius r as 

-„„,2 



(19) M{r) = 



1 — cos 77 



G 



(27) 
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4. Applications 

The procedure described above is applied to estimate a 
lower limit for the mass profile of the Coma cluster. For 
t his p urpose, we employ data provided by Geller et al. 
( |l999|) . The largest sample presented by these authors 
contains 691 galaxies. The cluster's velocity czc (c is 
the speed of light and zc is the redshift of the clus- 
ter centre) is 7090Kmsec^^ and for a Hubble constant 
Ho = 100hKmsec"iMpc~\ its distance is 70.9h-iMpc. 

As noted in Sect. 2, our approach is based on the pro- 
file of the amplitude A of the caustics. Departures from 
spherical symmetry, random motions due to the develop- 
ment of substructure, the finite number of galaxies in the 
cluster and observational errors are some of the reasons 
making the determination of caustics non-trivial. A gen- 
eral approach, given by Diaferio (1999), is based on the ar- 
gument that caustics are high density regions (see Regos & 
Geller 1989) and is as follows. First, a smooth estimate for 
the density / of observed galaxies on the (i?, Vobs) plane is 
calculated and then a cut is applied at some density con- 
tour which is taken to correspond to the caustics. A sim- 



ilar procedure is followed by Reisenegger et al. (2000) in 
their application to the Sharpley Supercluster. Problems 
associated with the above procedure, such as the choice 
of smoothing lengths and the density cutoff, are discussed 
in the above p apers . The dashed line -shown in Fig. ^ is 
Geller et al.'s ( 199E ) estimation of the amplitude Aest{R) 
of the caustics derived by the above-discussed procedure. 
We assume that this can be written in the form 

AestiR) ^ A{R) + dA{R) 

where A{R) stands for the component due to the radial 
infall velocity and the positive dA(R) describes the 
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Fig. 2. Dashed line: profile of the amplitude of the caus- 
tics (Geller et al. ( |1999| )). Solid curve: approximation of 
this line using analytical smooth functions (see eqn. 
and text for more details.) 
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Fig. 3. Radial velocity profile corresponding to the 
smooth curve plotted in Fig.|| 



'inflation' of the caustics due to the development of ran- 
dom motions. Unfortunately, dA(R) cannot be modelled 
properly and an infinite number of curves lower than A^st 
could be provided as models for A. For this reason we use 
the following procedure. First, we use models for A{R) of 
the form 



AiR) = vo 



R 



Wo > 0, R< ru- 



(28) 



Then, we estimate the values of the parameters vo,rM 
and k in order that these curves satisfy the condition 
dA{R) > for all values of R and minimize the sum 
"^[AestiR) — A{R)]'^. In this way, we derive a profile that 
could be considered as an upper bound for the A compo- 
nent. Such a curve is plotted in Fig.^. The parameters used 
in eqn.dl) are: vq = 2 TOOKmsec^i, ru = 60ft.-^Mpc, k = 
2.8. We note that ru is not the physical radius of the 
system, but a fitting parameter. The form of eqn. (28) 
has mathematical advantages, since it allows an analyt- 
ical evaluation of the profile of the infall velocity, using 
eqns. (^) and (||) reported in Sect. 2. This is given by the 
relation 



Vr{r) 



-Vo 



R. 



2 n 

2k-l 



r<rM, k>^ (29) 



with r = 

In Fig.0, the resulting profile of the radial infall veloc- 
ity is plotted. Assuming that the above profile is a good 
approximation of the radial systematic component of the 
infall velocity, we apply the spherical infall model to calcu- 
late mass profiles, for various values of the density param- 
eter ilo of the universe, that are plotted in Fig.^. Starting 
from the bottom of the figure, the different curves cor- 
respond to ilo — 0.2,0.4,0.6,0.8,1.0 respectively. Larger 
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Fig. 4. Mass profiles corresponding to the radial veloc- 
ity profile (Fig.^) for varying fJo. From the bottom of the 
figure, the curves correspond to flo — 0.2,0.4,0.6,0.8,1.0 



respectively. Dashed curve: Navarro et al. ( 1995 ) mass pro- 
file. 
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Fig. 5. Mass profiles for various values of as in Fig 
The bar at 0.5 shows the limits given by Hughes (19891 
while the bar at 2.5 shows the limits given by Briel et al. 
(1992). The outer bar at 10 shows the estimation given by 
Geller et al. (|l999|) 



values of ilp lead to steeper mass profiles at the outer re- 
gions. The variation of Qq alters the profile at the outer 
regions of the cluster and (as expected) larger values of 
the density parameter of the universe give rise to larger 
amount of mass inside a given radius. Notice that the mass 
scale s as flo^''^, roughly as expected (see Regos & Geller 
198E ). Thus, the matter inside 10/i~^Mpc varies for vari- 



ous values of the density parameter of the universe in the 
range 0.994 to 1.12 x 10^^ H-^Mq. The dashed curve, plot- 
ted also in this Fi g., is t he cumulative mass of an halo with 
a Navarro et al. ( 1995| ) profile. This is given by the rela- 
tion M(r) = 47r/9o?'^s{ln[l -|- (r/rs)] - (r/rs)/[l -|- (r/rs)]}, 
where po and are fitting parameters. The values of 
the parameters used are: po — 3.8 x 10^^/i^MqMpc~^, 
Ts = 0.2/i~^Mpc. It is clearly shown in Fig.^that there is a 
remarkable agreement of the form of the cumulative mass 
profiles estimated by t he pre sent approach (solid curves) 
with the Navarro et al. (1995) mass profile (dashed curve). 
Note that the dashed curve describes the mass profile in 
the virialized region of the cluster. 

5. Conclusions 

Several methods have been proposed in the literature re- 
garding the estimation of the mass of galaxy clusters. 

A class of such methods, used to estimate the den- 
sity profile of these systems, is based on the use of X- 



ray data (The & White 198S; Hughes 1989; Watt et al 



1992). It is assumed that the gas is described by the equa- 



tion of hydrostatic equilibrium. The latter is solved using 
the density and temperature profiles of the gas, that are 
measured using X-ray observations. This method is accu- 
rate in the region where the hydrostatic equilibrium holds. 



that is t he cen tral region of galaxy clusters. For example, 
Hughes ( 1989|) determined the minimum and maximum 
mass profiles consistent with all the data for Coma -if 
the mass-to-light ratio of the data was not constrained to 
be a constant- obtaining allowed mass ranges of M0.5 = 
1.9-3.6xlO"/i-iM0 andMa.s = 0.5- 1.5 x lO^^/i-iM© 
where M0.5 and M2.5 are the masses inside 0.5ft.~"'^Mpc 



and 2.5/i~^Mpc respectively. Briel et. al (1992) used their 
ROSAT survey image to measure the X-ray surface bright- 
ness out to 100' from the cluster center. They found that 
the binding mass is more centrally concentrated than the 
X-ray gas, and obtain M2.5 = 0.6 - 1.2 x W^hr^MQ. 

Another class of methods employs the projected veloc- 
ities of the galaxies-members of the cluster. The methods 
based on the form of the caustics in redshift space belong 
-among others- in this class. In a series of papers (Diaferio 
et al. [1997| Diaferio |1999| ; Geller et al. |1999| ), it has been 
proposed to estimate the mass profile of a galaxy cluster 
using the equation 



GM(t) « - / ^Lt(a;)da:. 



(30) 



Although there is no rigorous proof for this relation, its 
testing against cosniological N-body simulations shows 
that it approximates well the mass profiles. The deriva- 
tion of the above relation is based on the argument that 
the amplitude of the velocity field of galaxy clusters de- 
pends mainly on local dynamics, since the random motion 
is significant. The mass profile is, in any case, roughly 
proportional to the square of the amplitude of t he ca us- 
tics and since the amplitude used by Geller at al. (199£) is 
larger because it includes the effect of random motion, one 
expects their estimation to result in a larger mass for the 
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system. They estimate a mass inside a radius of lOh ^Mpc 
in the range (1.65 ± 0.41) x lO^^h-^MQ. 

In Fig.|, we compare the mass estimations provided by 
our method for the Coma cluster to the mass estimations 



given by Hughes ( |l98g| ) , Briel (|l99g) and Geller et al. 
( [19991 ). It is clear that the mass profile resulting from our 
application shows the minimum estimated values of mass 
at different radii and provides a reasonable estimation of 
the lower limit of the mass profile of the system. 

More detailed observations could lead to a more accu- 
rate definition of the caustics and improve the estimation 
of mass profiles . 

Acknowledgements. Thanks to the EMPIRIKION Foundation 
for its support. 
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